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I , Abstract 

. The original B. Pontecorvo idea of neutrino oscillations is dis- 

\ cussed. Neutrino mixing and basics of neutrino oscillations in vacuum 

CNl ' and in matter are considered. Recent evidences in favour of neutrino 

. oscillations, obtained in the Super-Kamiokande, SNO, KamLAND and 

\ other neutrino experiments are discussed. Neutrino oscillations in the 

<^ ' solar and atmospheric ranges of the neutrino mass-squared differences 

r"| ' are considered in the framework of the minimal scheme with the mix- 

Q-i! ing of three massive neutrinos. Results of the tritium /3-decay experi- 

QhI ments and experiments on the search for neutrinoless double /3-decay 

^ ■ are briefly discussed. 

_> : 

[x^ ' 1 Pontecorvo idea of neutrino oscillations 

- - ' B. Pontecorvo came to the idea of neutrino oscillations in 1957, soon after 

parity violation in /3-decay was discovered by Wu et. al. P and the two- 
component theory of massless neutrino was proposed by Landau |2] , Lee and 
Yang and Salam 0. 

For the first time B. Pontecorvo mentioned the possibility of neutrino 
antineutrino transitions in vacuum in his paper on muonium ^ antimuo- 
nium transition {fi~e~^ ^ fi'^e~). He believed in the analogy between the 
weak interaction of leptons and hadrons and looked in the lepton world for a 
phenomenon that would be analogous to K'^ ^ K'^ oscillations. In his work 
of 1957 ISj B. Pontecorvo wrote: 

^Dedicated to B. Pontecorvo in 90th anniversary oj his birth. 
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"If the two- component neutrino theory should turn out to be 
incorrect (which at present seems to be rather improbable) and 
if the conservation law of neutrino charge would not apply, then 
in principle neutrino antineutrino transitions could take place 
in vacuum." 

B. Pontecorvo wrote his first paper on neutrino oscillations later in 1958 [H] 
at the time when F. Reines and C. Cowan [7^ had just finished the famous 
reactor experiment, which led them to discover the electron antineutrino 
through the observation of the inverse /^-process 

z7e + p — > e"*" + . (1) 

At that time R. Davis was doing an experiment with reactor antineutrinos 
[Hj. R.Davis searched for production of ^^Ar in the process 

z7,+ 37ci^e-+ 37Ar, (2) 

which is allowed only if the lepton number is not conserved. A rumor reached 
B. Pontecorvo that R. Davis had seen some events (2). B. Pontecorvo, who 
had earlier been thinking about possible antineutrino ^ neutrino transitions, 
decided to study this possibility in details. 

"Recently the question was discussed [S] whether there exist other 
mixed neutral particles beside the mesons, i.e., particles that 
differ from the corresponding antiparticles, with the transitions 
between particle and antiparticle states not being strictly for- 
bidden. It was noted that the neutrino might be such a mixed 
particle, and consequently there exists the possibility of real neu- 
trino ^ antineutrino transitions in vacuum, provided that lepton 
(neutrino) charge is not conserved. This means that the neutrino 
and antineutrino are mixed particles, i.e., a symmetric and anti- 
symmetric combination of two truly neutral Majorana particles 
Ui and 1^2 of different combined parity" 

In Ref. B. Pontecorvo wrote that he has considered this possibility 

". . . since it leads to consequences which, in principle, can be 
tested experimentally. Thus, for example, a beam of neutral lep- 
tons consisting mainly of antineutrinos when emitted from a nu- 
clear reactor, will consist at some distance R from the reactor of 
half neutrinos and half antineutrinos." 
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In 1958 only the electron neutrino was known. In Ref. [H] B. Pon- 
tecorvo considered oscillations of the active right-handed antineutrino to a 
right-handed neutrino, the only possible oscillations in the case of one type of 
neutrino. In this paper, that was written at a time when the two-component 
neutrino theory had just appeared and the Davis reactor experiment was not 
yet finished, B. Pontecorvo discussed a possible explanation of the Davis 
"events". Later B. Pontecorvo understood that in the framework of the two- 
component neutrino theory, established after M. Goldhaber et al. experi- 
ment [D], right-handed neutrinos are practically sterile particles.^ B. Pon- 
tecorvo was in fact the first, who introduced the notion of sterile neutrinos 
in Ref. [10]. 

B. Pontecorvo stressed [Hj that if the oscillation length R is large, it 
will be difficult to observe the effect of the neutrino oscillations in reactor 
experiments. He wrote: 

"Effects of transformation of neutrino into antineutrino and vice 
versa may be unobservable in the laboratory because of the large 
values of R, but will certainly occur, at least, on an astronomical 
scale." 

B. Pontecorvo came back to the consideration of neutrino oscillations 
in 1967 [lOj. At that time the phenomenological V-A theory of Feynman 
and Gell-Mann ^2), Marshak and Sudarshan ^2] was established and in 
the Brookhaven experiment which was proposed by B. Pontecorvo in 

1959 [TF, it was proved that (at least) two types on neutrinos z/g and i/^ exist 
in nature. B. Pontecorvo confidence in nonzero neutrino masses and neutrino 
oscillations became stronger after these findings. 

In Ref. PU] B. Pontecorvo discussed experiments in which conservation of 
lepton numbers were tested and came to the conclusion that there is "plenty 
of room for violation of lepton numbers, especially in pure leptonic interac- 
tions". He formulated jHlE] the conditions at which neutrino oscillations 
are possible: 

1. Lepton numbers, conserved by the usual weak interaction, are violated 
by an additional interaction between neutrinos. 

^If the lepton number is violated and neutrinos with definite masses are Majorana 
particles the process (j2J) in principle is allowed. However, the cross section of this process 
is strongly suppressed by the factor (m/E)'^ ( m is the neutrino mass and E the neutrino 
energy) 
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2. Neutrino masses are different from zero. 

In modern terminology these conditions are equivalent to the assumption 
that in the total Lagrangian there is a neutrino mass term, non-diagonal in 
flavour neutrino fields. 

In Ref. [TU] B. Pontecorvo considered different lepton numbers in the 
case of two types of neutrinos: additive electron and muon lepton numbers 
Le and L^, ZKM lepton number jTH] equal to 1 for e~ and and -1 for e+ 
and fi~ etc. He stressed that if violation of and occurs, in addition to 
the oscillations between active flavour neutrinos z/^ ^ oscillations between 
active and sterile neutrinos ( ^ u^l, ^ i^^iL etc) become possible. 

In the case of ZKM lepton number L for one four-component neutrino u 
we have the following correspondence: ul = Ue, J^r = i^^i, i^L = ^^i, VR = Ve- 
If the interaction between neutrinos does not conserve L, only transitions 
between active neutrinos {vl ^ vl and ^ vr) are then possible. 

In his work of 1967 [TO] B. Pontecorvo discussed the oscillations of so- 
lar neutrinos. At that time R. Davis started his famous experiment on the 
detection of the solar neutrinos. The experiment was based on the radio- 
chemical method, proposed by B. Pontecorvo in 1946 jTTj. Solar neutrinos 
were detected via the observation of the Pontecorvo-Davis reaction 

z/e + ^^Cl ^ e" + 3^Ar. (3) 

Quite remarkably, in llOjj B. Pontecorvo envisaged the solar neutrino prob- 
lem. Indeed, before the flrst results of Davis experiment were obtained, he 
wrote: 

"From an observational point of view the ideal object is the sun. 
If the oscillation length is smaller than the radius of the sun region 
effectively producing neutrinos, (let us say one tenth of the sun 
radius Rq or 0.1 million km for neutrinos, which will give the 
main contribution in the experiments being planned now), direct 
oscillations will be smeared out and unobservable. The only effect 
on the earth's surface would be that the flux of observable sun 
neutrinos must be two times smaller than the total (active and 
sterile) neutrino flux." 

It was shown by Gribov and Pontecorvo in 1969 ^7] that in the neutrino 
mass term only active left-handed flelds z/gL and v^i can enter. The corre- 
sponding mass term is called Majorana mass term. Particles with deflnite 
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masses are in this case Majorana neutrinos and flavour fields Uei and u^l are 
connected with the Majorana fields and V2L by the mixing relation 

Pe = cos 6' viL + sin 6* 1^21] 

(4) 

z/^ = - sin 9 + cos 9 P2L, 

where 9 is the mixing angle. 

In Ref. jT7j the expression for the survival probability of z/g was obtained 
and oscillations of solar neutrinos in vacuum were considered. 

In the seventies the Cabibbo-GIM mixing of d and s quarks was estab- 
lished. In Ref. [T^ neutrino oscillations in the scheme of the mixing of two 
Dirac neutrinos ui and 1^2 were considered. This scheme was based on the 
analogy between quarks and leptons. In the same work possible values of 
the mixing angle 9 were discussed: 

". . . it seems to us that the special values of mixing angle 9 = 
(the usual scheme in which muonic charge is strictly conserved) 
and 9 = tt/4 are of the greatest interest^. " 

Neutrino oscillations and in particular solar neutrino oscillations were 
considered in Ref. in the general case of the Dirac and Majorana mass 
term: the latter does not conserve any lepton numbers. It is built from 
active left-handed fields uil and sterile right-handed fields uui. 

In 1978 the first review on neutrino oscillations was published ^20)j. At 
that time the list of papers on neutrino oscillations was very short: except 
for the cited papers of B. Pontecorvo and his collaborators, it included only 
references [2111^1231211123 and j2nilZZ|- The review attracted the attention 
of many physicists to the problem of neutrino mass and neutrino oscillations.^ 
In particular, it was stressed in this review that due to the interference 
character of neutrino oscillations the investigation of this phenomenon is the 
most sensitive method to search for small neutrino mass-squared differences. 

During the preparation of the review [201 the authors became familiar 
with the work of Maki et al. [21] , which was unknown to many physicists 

■^We know from today's data that small and 7r/4 values of neutrino mixing angles are 
really the preferable ones. Indeed, as we will see later, in the framework of the mixing of 
three massive neutrinos with three mixing angles the angle 612 is small and the angles ^13 
and 023 are close to 7r/4. 

^It has now about 500 citations. 
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for many years. In this paper, in the framework of the Nagoya model, in 
which p, n and A were considered as bound states of some vector particle 

and leptons, the mixing of two neutrinos was introduced. In addition to 
the usual neutrinos Ve and z/^, which they called weak neutrinos, the authors 
introduced massive neutrinos vi and V2i which they called true neutrinos. In 
the Nagoya model the proton was considered as a bound state of the 
and the neutrino Vi. 

Maki, Nakagawa and Sakata assumed that the fields of weak neutrinos 
and of true neutrinos are connected by an orthogonal transformation: 

"It should be stressed at this stage that the definition of the par- 
ticle state of the neutrino is quite arbitrary; we can speak of neu- 
trinos which are different from weak neutrinos but are expressed 
by a linear combination of the latter. We assume that there exists 
a representation which defines the true neutrinos through some 
orthogonal transformation applied to the representation of weak 
neutrinos: 

vi = +Uf, cos 5 + sin 5 

(5) 

1^2 = —Ve sin 5 + COS 5 " 

In Ref. [21] neutrino oscillations, as a phenomenon based on the quantum 
mechanics of a mixed system, were not considered. In connection with the 
Brookhaven experiment ^3], which was going on at that time, the authors 
wrote 

"In the present case, however, weak neutrinos are not stable due 
to occurrence of virtual transmutation v^, induced by the 

interaction C = gV2h'2X^X (X is a field of heavy particles). If 
the mass difference between ui and z/2, i.e. \m^^ — my^\ = m^^, 
is assumed to be a few MeV the transmutation time T {Uf. ^ v^) 
becomes ~ 10^^^ sec for fast neutrinos with momentum ~ BeV/c. 

Therefore a chain of reactions such as 

1^^ + Z {nucleus) Z' + (^u^ / e~) (6) 
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is useful to check the two-neutrino hypothesis only when \m^^ — 
i^uil ^ 10^^ MeV under the conventional geometry of the exper- 
iments. Conversely, the absence of e~ in the reaction ® will be 
able not only to verify two-neutrino hypothesis but also to pro- 
vide an upper limit of the mass of the second neutrino 1^2 if the 
present scheme should be accepted." 

Recently very strong, convincing evidence in favour of neutrino oscilla- 
tions, which we will discuss later, were obtained. It required many years of 
work and heroic efforts of many experimental groups to reveal the effects 
of tiny neutrino masses and neutrino mixing. From our point of view it is a 
proper time to give a tribute to the great intuition of B. Pontecorvo, who 
pursued the idea of neutrino oscillations for many years at a time when the 
general opinion, mainly based on the success of the two-component neutrino 
theory, favoured massless neutrinos ^. The history of neutrino oscillations is 
an illustration of the importance of analogy in physics. It is also an illus- 
tration of the importance of new courageous ideas, not always in agreement 
with the general opinion. 

2 Introduction 

Convincing evidence of neutrino oscillations was obtained in the last 4- 
5 years in experiments with neutrinos from natural sources: in the atmo- 
spheric 123 EOl EH and solar neutrino experiments |H21 El ESI EH EH] • 
Recently strong evidence in favour of neutrino oscillations was obtained in 
the reactor KamLAND experiment j39j . 

The observation of neutrino oscillations implies that neutrino masses are 
different from zero and that the states of flavour neutrinos z/g, i/^, Uj. are 
coherent superpositions (mixtures) of the states of neutrinos with definite 
masses. All existing data, including astrophysical ones, suggest that neu- 
trino masses are much smaller than the masses of leptons and quarks. The 
smallness of neutrino masses is a signature of a new, beyond the Standard 
Model, physics. 

In this review we present the phenomenological theory of neutrino masses 
and mixing and the see-saw mechanism of neutrino mass generation jlUj. We 
consider neutrino oscillations in vacuum and transitions between different 

^For a collection of the papers of B. Pontecorvo see Ref. |28|. 
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flavour neutrinos in matter. Then we will discuss the recent experimental ev- 
idence in favour of neutrino oscillations: the results of the Super-Kamiokande 
atmospheric neutrino experiment in which a significant up-down asym- 
metry of the high-energy muon events was observed, the results of the SNO 
solar neutrino experiment (3^1 EZl IHH] ? in which a direct evidence for the 
transition of the solar z/g into z/^ and Uj. was obtained and the results of 
the KamLAND experiment [21] in which the disappearance of the reactor 
Ue was observed. We will also discuss the long baseline CHOOZ jHj and 
Palo Verde jj^ reactor experiments in which no indications in favour of neu- 
trino oscillations were obtained. Yet the results of these experiments are very 
important for neutrino mixing. 

Neutrino oscillations in the atmospheric and solar ranges of neutrino 
mass-squared differences will be considered in the framework of three-neutrino 
mixing. We will show that neutrino oscillations in these two regions are 
practically decoupled and are described in the leading approximation by the 
two-neutrino formulas. 

The investigation of neutrino oscillations is based on the following as- 
sumptions, which are supported by all existing experimental data. 

1. The interaction of neutrinos with other particles is described by the 
Standard Model of the electroweak interaction. 

The Standard Charged Current (CC) and Neutral Current (NC) La- 
grangians are given by 

Here g is the SU(2) gauge coupling constant, 6w is the weak angle, 
and are the fields of the charged (W^^) and neutral (Z°) vec- 
tor bosons and the leptonic charged current j^^ and neutrino neutral 
current j^'^ are given by the expressions: 

2. Three flavour neutrinos v^, and exist in nature. 

From the LEP experiments on the measurement of the width of the 



decay Z v ^ v the following value was obtained for the number 
of flavour neutrinos n,. 
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n^^ = 3.00 ± 0.06 (9) 
and from the global fit of the LEP data the value 

n^^ = 2.984 ± 0.008. (10) 

was found. 

3 Neutrino masses and mixing 

If neutrino fields enter only in the Standard Model (SM) Lagrangians (|7j), 
neutrinos are massless particles and the three flavour lepton numbers, elec- 
tron (Lg), muon (L^) and tau (Lt), are separately conserved. The hypothesis 
of neutrino mixing is based on the assumption that in the total Lagrangian 
there is a neutrino mass term, which does not conserve flavour lepton num- 
bers. Several mechanisms of the generation of the neutrino mass term were 
proposed. In this section we will discuss the most plausible see-saw mecha- 
nism ||40j. 

Let us start with the phenomenological theory of the neutrino masses and 
mixing. There are two types of possible neutrino mass terms (see jiil liH]). 

1. Dirac mass term 

= -jy'j^M^iy'^ + h.c. (11) 

Here 



( VeL \ 



V 



R 



( J^eR \ 
I^^lR 
VrR 



(12) 



and M° is a complex non-diagonal matrix. We assume that in the 
column z/^ not only left-handed flavour neutrino fields z/gL, v^l^ ^tL but 
also sterile fields could enter. Sterile fields do not enter into the stan- 
dard charged and neutral currents (see discussion later). 
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The matrix M can be diagonalized by a bi-unitary transformation. 
We have (see 

M° = 1/mf/+, (13) 

where V and U are unitary matrices and rriik = rriiSik', rrii > 0. From 
Eq. (jllj) and Eq. (fT^ for the neutrino mass term we obtain the standard 
expression 

= ^J^m^UiUi, (14) 

i 

where Ui is the field of neutrino with mass mj. 

The flavour field uil is connected with left-handed field by the 
mixing relation 

^lL = ^UuViL. (15) 

i 

It is obvious that the total Lagrangian is invariant under the global 
gauge transformation 

where a is an arbitrary constant phase. From this invariance it follows 
that the total lepton number L = Li is conserved and z/j is the field 
of the Dirac neutrinos and antineutrinos with L = 1 for neutrinos and 
L = —1 for antineutrinos. 

2. Majorana mass term 

£^ = -i(^)^M^z.i + h.c. (16) 

Here 

where C is the charge conjugation matrix, which satisfies the conditions 

The matrix is a symmetrical matrix. In fact, taking into account 
the Fermi-Dirac statistics of the field z/^, we have 
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From this relation we obtain 



M 



■M 



For a symmetrical matrix we have (see |m I45j) 

= {U+fmU+, 



(17) 



where U is a unitary matrix and rriik = rriidik', rrii > 0. 

From Eqs.()16|) and ()17p it follows that the Majorana mass term takes 
the form 



where z/j is the field of neutrino with mass mj, which satisfies the Ma- 
jorana condition 



The flavour field un is connected with the Majorana fields by the 
mixing relation 



In the case of the Majorana mass term there is no global gauge in- 
variance of the total Lagrangian. Hence, Majorana neutrinos are truly 
neutral particles: they do not carry neither electric charge nor lepton 
numbers. In other words Majorana neutrinos and antineutrinos are 
identical particles. 

If there are only flavour fields un in the column i/^, the number of the 
massive neutrinos z/j is equal to the number of flavour neutrinos (three) and 
?7 is a 3x3 Pontecorvo-Maki-Nakagawa-Sakata (PMNS) GDI IHl unitary 
mixing matrix. 

If, in the column i^^, there are also sterile fields Usl, the number of massive 
neutrinos i^i will be larger than three. In this case the mixing relation takes 
the form 






(19) 




(20) 
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3+ns 3+ns 

I^IL = ^ UiiUiL] VsL = Usil^iL , (21) 

i=l i=l 

where is the number of the sterile fields, f/ is a unitary (3 + n^) x (3 + n^) 
matrix, i/j is the field of neutrino with mass rrii {i = 1, 2, .., 3 + n^). 

Sterile fields can be right-handed neutrino fields, SUSY fields etc. If more 
than three neutrino masses are small, the transition of the flavour neutrinos 
^e, ^fi, into sterile states becomes possible. 

A special interest has the so-called Dirac and Majorana mass term. Let 
us assume that, in addition to the flavour fields un, in the column z/^ also 
the fields {i^irY {I = e,/i, r) enter. ^ The Majorana mass term (fT^ can then 
be presented in the form of the sum of the left-handed Majorana, Dirac and 
right-handed Majorana mass terms: 



Here 



= -IWM^^^'u'. + h.c. (22) 



(23) 



Mj^ and are complex non-diagonal symmetrical 3x3 Majorana matrices 
and is a complex non-diagonal 3x3 Dirac matrix. After the diagonal- 
ization of the mass term ()22|) one gets 









( VeR \ 
















^ VrR / 



VlL = UhU.l; [virY = UuViL, (24) 
i=l i=l 

where U is the unitary 6x6 mixing matrix and fi is the field of the Majorana 
neutrino with mass m,-. 
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(vrY is the left-handed component of the conjugated field v'^. In fact, we have 

{^hY-^-^cd^^uI. 
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The standard see-saw mechanism of the neutrino mass generation jlHI is 
based on the assumption that the neutrino mass term is the Dirac and Ma- 
jorana one, with = 0. To illustrate this mechanism we will consider the 
simplest case of one type of neutrino. Let us assume that the standard Higgs 
mechanism with one Higgs doublet, which is the mechanism of generation 
of the masses of quarks and leptons, also generates the Dirac neutrino mass 
term 

£D = -mz/Ri/L + h.c. (25) 

It is natural to expect that the mass m is of the same order of magnitude 
of the lepton or quark masses in the same family. However, we know from 
experimental data that neutrino masses are much smaller than the masses 
of leptons and quarks. In order to "suppress" the neutrino mass we will 
assume that there is a beyond the SM, lepton number violating mechanism 
of generation of the right-handed Majorana mass term 

£^ = -iMz7«M^ + h.c., (26) 

with M > m (usually it is assumed that M ~ Mqut - lO^^GeV). The 
total mass term is then the Dirac and Majorana one, with 

- ( : s ) - ( ) 

After the diagonalization of the mass term one gets 



ul = i cos 6 uiL + sin 9 V2L 
(i^rY = -i sm6 + cos6 iy2L, (28) 

where ui and z/2 are the fields of the Majorana particles, with masses 

1 1 ^2 



mi = — M + - + 4m2 ^ — < m; 
1^. 1^ 



2.2 M ' (29) 



rras = - M + - + 4 ~ M 



The mixing angle 6 is given by the relation 



2 TTl 

tan2^ = — <1. (30) 
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Thus, the see-saw mechanism is based on the assumption that, in addi- 
tion to the standard Higgs mechanism of generation of the Dirac mass term, 
there exists a beyond the SM mechanism of generation of the right-handed 
Majorana mass term, which changes the lepton number by two and is char- 
acterized by a mass M ^ mJ The Dirac mass term mixes the left-handed 
field ul, the component of a doublet, with a singlet field {i'rY- As a result 
of this mixing the neutrino acquires Majorana mass, which is much smaller 
than the masses of leptons or quarks. 

In the case of three generations, in the mass spectrum there will be 
three light Majorana masses, much smaller than the masses of quarks and 
leptons, and three very heavy Majorana masses of the order of magnitude of 
the lepton number violation scale M (see ^OJISIEH])- 

Let us stress that if neutrino masses are of see-saw origin then: 

• Neutrinos with definite masses are Majorana particles. 

• There are three light neutrinos. 

• Heavy Majorana particles must exist. 

The existence of the heavy Majorana particles, see-saw partners of neu- 
trinos, could be a source of the barion asymmetry of the Universe (see |49j). 

4 Neutrino oscillations in vacuum 

In this section we will discuss the phenomenon of neutrino oscillations (see, 
for example, jSl ^1)- If in the total Lagrangian there is a neutrino mass 
term the flavour lepton numbers Le, L^, Lr are not conserved and transitions 
between different flavour neutrinos become possible. 

Let us first define the states of flavour neutrinos z/g, z/^ and u^- in the case 
of neutrino mixing. The flavour neutrinos are particles which take part in 
the standard weak processes. For example, the neutrino that is produced 
together with /i+ in the decay tt^ + is the muon neutrino z/^; the 

electron antineutrino z/g produces e+ in the process z^e + P — > e+ + n, etc. 

In order to determine the states of the flavour neutrinos let us consider a 
decay 

a-^h + l+ + (31) 
^It is obvious that for charged particles such mechanism does not exist. 
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If there is neutrino mixing, the state of the final particles is given by 

|/>= |/+) {ith\S\a), (32) 

i 

where is the state of neutrino with momentum p and energy 



and b \S\ a) is the element of the S-matrix. 

We will assume that the neutrino mass-squared differences are so small 
that the emission of neutrinos with different masses can not be resolved in 
the neutrino production (and detection) experiments. In this case we have 

(«/+6|5|a)~f/;* {uil+b\S\a)sM, (33) 

where {uil^blSl a)sM is the SM matrix element of the process (jSH), calcu- 
lated under the assumption that all neutrinos Vi have the same mass. From 
()32|) and ()33|) we obtain the following expression for the normalized state of 
the flavour neutrino vi: 

\^i) = Y.^u (34) 

i 

Thus, in the case of mixing of neutrino fields with small mass-squared differ- 
ences, the state of a flavour neutrino is a coherent superposition (mixture) of 
the states of neutrinos with definite masses.* 

In the general case of active and sterile neutrinos we have 

W) = "^Ki Wi) , (35) 

i 

where the index a takes the values e, /x, r, si, .... From the unitarity of the 
mixing matrix it follows that 

{j^a'Wa) = Sa'a- (36) 

The phenomenon of neutrino oscillations is based on the relation ()35|) . 
Let us consider the evolution of the mixed neutrino states in the vacuum. 

^The relation (|34|1 is analogous to the relation which connects the states of K'^ and 
K'^ mesons, particles with definite strangeness, with the states of Kg and mesons, 
particles with definite masses and widths. 
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If at the initial time t = a flavour neutrino z/^ is produced, the neutrino 
state at a time t will be: 



(37) 



where Ho is the free Hamiltonian. 

It is important that the phase factors in ()37p for different mass compo- 
nents are different. Thus, the flavour content of the final state jt'a)* will differ 
from the initial one. As we shall see later, in spite of the small neutrino 
mass-squared differences, the effect of the transition of the initial neutrino 
into another flavour (or sterile) neutrino can be large. 

Neutrinos are detected through the observation of CC and NC weak pro- 
cesses. By developing the state \i'a)t over the total system of neutrino states 
li'a), we have 



where 



-iHot 



,/,; e 



-iEA 



u*. 



(38) 
(39) 



is the amphtude of the transition z/q, z/q,/ during the time t. 

Taking into account the unitarity of the mixing matrix, from (|39p we 
obtain the following expression for the probability of the transition 
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Z/rv' 



Pfz/. 



Z^rv', 



5a'a + ^ Ua'iU*i (c" 



"'"jl 2E 



1) 



i>2 



(40) 



where L ~ t is the distance between a neutrino source and a neutrino detec- 
tor, E is the neutrino energy and Am^^ = mf — m\. 

Analogously, for the probability of the transition z/q, — ^ Va' we have: 



Pfz7„ 



Z^rv', 



(41) 



Let us stress the following general properties of the transition probabili- 



ties: 



^We label neutrino masses in such a way that mi < m2 < < 
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• Transition probabilities depend on L/E. 

• Neutrino oscillations can be observed if the condition 

is satisfied for at least one value of i. 

• From the comparison of Eqs. ()40|) and ()41|) we conclude that the fol- 
lowing relation holds 

This relation is the consequence of the CPT invariance, intrinsic for 
any local field theory. 

• If CP invariance in the lepton sector holds, the mixing matrix U is real 
in the Dirac case. In the Majorana case the mixing matrix satisfies the 
condition 

Uai = rj. , (42) 

where 77, = ±i is the CP parity of the Majorana neutrino z/j. From 
Eqs. (HO)), PTj) and we can conclude that in the case of CP invari- 
ance in the lepton sector we have the following relation 



In conclusion let us introduce the parameters which characterize the 3 x 
3 neutrino mixing matrix. In the general n x n case the unitary matrix 
is characterized by n{n — l)/2 angles and n{n + l)/2 phases. The phases 
of lepton fields are arbitrary. If neutrinos with definite masses are Dirac 
particles the phases of neutrino fields Ui are arbitrary as well. In this case 
the mixing matrix is characterized by 

n(n + l) . . . , , (n-l)(n-2) 
2[n — 1) — 1 (common phase) = 



physical phases. 

Thus, the 3x3 Dirac mixing matrix is characterised by three angles and 
one phase. Let us parameterize such matrix. Taking into account the unitary 
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condition |f/eiP = 1, we can choose the first two elements of the first row 
in following way 

Uel = l^esP COS ^12; f/e2 = V^l - [f/esP sin ^12 , (43) 

where 612 is the mixing angle. The third element of the first row can be 
parameterized as follows 

Ue3 = smen e-'\ (44) 

where 6 is the CP-phase and ^^13 the second mixing angle. Thus, for |i^e) we 
have 

We) = Vl- If/esP M + f/e3 Ws) , (45) 

where 

lz/12) = cos ^12 + sin ^12 1^2) ■ (46) 

The elements of the second and third line of the mixing matrix must be 
chosen in such a way that the condition 

{i^i'Wi) = 5n. (47) 

is satisfied. 

It is obvious that the two orthogonal vectors 

I'^i2) = ~ sin^i2 |z^i) + cos^i2 |z^2), (48) 

We) = -Ut, ki2) + ^l-\Ue,? ka) , (49) 

are orthogonal to the vector We)- 

We have introduced two angles ^12, ^13 and the phase 5. The third mixing 
angle ^23 will enter into play as follows 

I z/^) = cos 6123 |z/i^) + sin6'23 Wi) , (50) 

\v^) = - sin 6*23 + cos 6123 Wi) ■ (51) 

From Eqs. (jl^jl . (jMIj) and (j3T|) we can obtain all the mixing matrix elements. 
We shall be interested in the elements of the first line and third column. The 
former are given by Eqs. (j43|) and pijl . For the elements and Urs of 
the third column we have 
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f/^S = Vl - If/esP Sin^ss; f/r3 = a/1 " |f/e3p 008^23 (52) 

If the Vi are Majorana fields, their phases are fixed by the Majorana 
condition, Eq. p9|) . The Majorana mixing matrix can be presented in the 
form 

= t/^ (53) 

where is the Dirac mixing matrix and Sik [a) = e^("'-°i) 5ik is the phase 
matrix. 

The number of physical phases in the Majorana mixing matrix is equal 

to 

- -^K{n-1) = - (54) 

Thus, the 3x3 unitary Majorana mixing matrix is characterized by three 
mixing angles and three CP phases. 

From Eqs. ()39|) and ()53p it is easy to derive the following relation be- 
tween the amplitudes of the transition Ua — > I'a' for the Majorana and Dirac 
neutrinos: 



(55) 

Hence, the investigation of neutrino oscillations does not allow one to dis- 
tinguish between the cases of Dirac or Majorana neutrinos [HI]. 



5 Oscillations between two types of neutrinos 
in vacuum 

We will consider here the simplest case of the transition between two types 
of neutrinos. In this case the index i in Eq. ()40j) takes only one value {i = 2) 
and for the transition probability we have 

Pil^a ^ l^a') = \Sa'a + U^,2K^ (g-^^^'^l: - 1)|2 , (56) 

where Am? = m\ — m\. 
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From (j56p we obtain the following expression for the appearance proba- 
bility: 

P(z/a ^ z/q,) = -A„/;q (1 - cos Am^— ) (a' 7^ a), (57) 

where 

Aa';Q, = 4 |f/a'2p \Ua2\^ = P^a;a' ■ (58) 

Let us introduce the mixing angle 9. In the 2x2 case the mixing matrix 
U can be chosen to be real. Hence we set: 

Ua,2 = Sm6; Ua'2 = COS 6*. 

The amplitude Aa'-^a is then given by 

Aq';q = sin^ 29 . 

and the two-neutrino transition probability takes the standard form 

P(z/„ ^ z/„0 = - sin2 2^ (1 -cosAm^-^); (a' ^ ") (59) 
2 2E 

From Eqs. (j57|) and (|58|) it is obvious that in the two-neutrino case the 
following relations are valid: 

P(z/q -> Ua') = P(z/a' -> l^a) = P{i^a i^a') ■ (60) 

Thus, the CP violation in the lepton sector can not be revealed in the case 
of transitions between two types of neutrinos. 

The survival probability P(z/q, z/q,) is determined by the condition of 
conservation of probability: 

P(z/„ ^ i/„) = 1 - P(z/« ^jy^,) = 1-- sin^ 26 (1 - cos Am^— ) . (61) 

2 2E 

From Eqs. (jUUI) and (pT|) it follows that the two-neutrino survival prob- 
abilities satisfy the following relation: 

P{Ua fa) = P{l^a' ^ l^a') ■ (62) 

Thus, only two oscillation parameters, sin^ 26' and Am^, characterize all 
transition probabilities in the case of a transition between two types of neu- 
trinos. 
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The expressions and (jHH) describe periodic transitions between two 
types of neutrinos (neutrino oscillations). They are widely used in the anal- 
ysis of experimental data.^° 

Let us notice that the expression for the two-neutrino transition 
probability can be recast in the form 

P(z/^ i^a') = - sm^29 (1 - cos2 7r^) , (63) 
2 Lq 

where 

Lo = 4 7r-^ (64) 

is the oscillation length. 

Finally, the two-neutrino transition probability and the oscillation length 
can be rewritten as 

P(z/„ i^„/) = - sin^ 29 (1 - cos 2.53 Am^^) (^5) 
2 E 

and 

where E is the neutrino energy in MeV (GeV), L is the distance in m (km) 
and Am? is the neutrino mass-squared difference in eV^. 



6 Neutrino oscillations in matter 

In the previous sections we have considered neutrino oscillations in the 
vacuum. However, if neutrinos pass through the Sun, the Earth, supernova, 
etc., matter can significantly alter the neutrino mixing and the probabilities 
of the transitions between different types of neutrinos. Here we will discuss 
the effects of matter on the neutrino transition probabilities (see |52| 13^ 13^ 

113 ESI). 

The refraction index of a particle with momentum p is given by the fol- 
lowing classical expression 

^"^As we shall see later, in the case of three-neutrino mixing neutrino oscillations in 
different ranges of Am^ are described, in leading approximation, by the two-neutrino 
formulas. 
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^^=1 + ^/(0)/^' (67) 

where p is the number density of matter and /(O) is the amphtude for 
elastic scattering in the forward direction. The second term in Eq. ()65|) is 
due to the coherent scattering of the particle in matter. 

In the case of neutrinos the amplitude of the process z/gC — > i/gC is different 
from the amplitude of the processes f^^r^ t'/i.re. This is related to the 
fact that to the matrix element of the process z/gC — > z/gC give contribution 
diagrams with the exchange of the W and Z bosons while to the matrix 
element of the processes v^^t^ i^^j.,t^ only the diagram with the exchange 
of the Z boson contributes. Thus, the refraction indexes of and v^^^ are 
different. If neutrino mixing occurs, this difference between the refraction 
indexes leads to important effects for the neutrino transitions in matter. 

In the flavour representation the evolution equation of a neutrino with 
momentum p in matter has the following general form: 

* = E(^o + ^^)-' (68) 

a' 

Here Hq is the free Hamiltonian, Hi is the effective Hamiltonian of the inter- 
action of neutrino with matter and ^^(t) is the amplitude of the probability 
to find Va at a time t. We have 



{H^)aa' =<«|Ho|a'>= J]U,iEiU:,i, (69) 



where U is the neutrino mixing matrix in vacuum and E'j ~ p + 

Let us consider the effective Hamiltonian of the interaction of neutrinos 
with matter. We will limit ourselves to the case of the flavour neutrinos. 

Due to the v^ — v^j^ — Vr universality of the NC, the NC part of the effective 
Hamiltonian is proportional to the unit matrix. This part of the Hamiltonian 
can be excluded from the evolution equation through a phase transformation 
of the function a(t). Thus, we need to take into account only the CC part of 
the Hamiltonian of z/g — e interaction. We have: 

{lii)ee =< z/emat| ^ 2 z7eL7"z/eL 67^(1 - 75)e |z/emat >~ V2GFPe , (70) 
V2 
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where Iz/ginat > is the state vector of a neutrino with momentum p and 
matter and is the electron number density. 

We will consider now the simplest case of two flavour neutrinos. The 
mixing matrix has the form 

/ cos6' sin^ \ , 
\ — sm U cos 6 1 

where 9 is the vacuum mixing angle. From Eqs. (jU^ -()7H ) the total effective 

Hamiltonian of neutrinos with momentum p in matter can then be written 
as:^2 

1 / -Am^cos26' + v4 Am^sin26'\ , . 

4^1^ Am2sin20 Am^ cos 2^- A J ^ ' 

where 

Am^ = m\ — m\] A = 2\/2GFPeP- 
The total Hamiltonian ()72|) can be easily diagonalized. We have 

^ ^jjm^m jjm+ ^ ^73^ 

where E"^ = 1, 2) are the eigenvalues of the matrix H and 

/ cos^'" sin^-X 

V -sin^*" cosO"' ) ^ ' 

is the orthogonal mixing matrix in matter. 

For the eigenvalues we find the following expressions: 

= ^l.^(Am2 cos 20 - Af + {Atv? sin 29^), (75) 
while the mixing angle 6"^ is given by the relations 
tan 2^- Am^sin2^ 



Am2 cos 26- A' 

Am^ cos 26 - A , 
cos 26"' = , . 76 

v/(Am2 cos 26 - Ay + (Am^ sin 26^ 



^^The same expression for the effective Hamiltonian of neutrinos in matter can be ob- 
tained from Eq. (see [S]) 

^^We omit the irrelevant unit matrix. 
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Let us consider now the simplest case of constant electron density. The 
probability of the transition ^a' 7^ «') in niatter is given in this case 

by the expression 

P'"(z/« ^ = P'"(z/„, ^ z/^) = i sin^ 20™ (1 - cos AE™ L) , (77) 
where L is the distance which the neutrino travels in matter and 

AE™ = E™ - = — v/(Am2 cos 20 - A) 2 + (Am2 sin 26 f . (78) 

2p 

From Eqs. ()77p and ()78|) it follows that the oscillation length in matter is 
given by 

° v/(Am2 cos 20 - AY + (Am2 sin 20)2 ' 

By putting pe = 0, we can easily see that the expressions (f77j) and (f?^ 
reduce to the corresponding vacuum expressions (j63|l and (jHll)- Again the 
survival probability of z/q {va') is given by the condition of conservation of 
the total probability: 

P'^il^a ^ l^a) = P'^i^a' ^ ^a') = 1 " P'^i^a ^ ^a') ■ (80) 

The neutrino mixing angle and oscillation length in matter can be signif- 
icantly different from the vacuum values. Let us assume that the condition 

Am^cos 20 = A = 2 v^G^PeP (81) 

is satisfied. Then from ()76p it follows that the mixing in matter is maximal 
(0™ = 7r/4), independently on the value of the vacuum mixing angle 0. In 
this case the oscillation length in matter turns out to be: 



L 



m 



'° sin 20 
where Lq is the vacuum oscillation length. 

In general the electron density is not a constant. For example, in the Sun 
the density is maximal in the center of the Sun and it decreases practically 
exponentially to its periphery. 

Hence in the general case the evolution equation of neutrinos in matter 
has the form 

.da{t) 



dt 
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i7™(t) ait) , (82) 



where H"^{t) = Hq + Hj{t) is the total effective Hamiltonian in the ffavour 
representation. 

This equation can be easily solved in the case of slowly changing density 
of electrons pe- In fact, the hermitian Hamiltonian H'^{t) can be diagonalized 
by the unitary transformation 

H^'lt) = U'^it) E'^it) W^+it), (83) 

where U^{t)U'^^{t) = 1 and E^{t) is the eigenvalue of the Hamiltonian. Let 
us now introduce the function 

a'(t) = t/™+(t)a(t); (84) 

from Eqs. (IH2I)-(|H1|) we have 

^ ^ = (E-^it) - ^ U-\t) ^^) a\t) . (85) 

If we assume that the function P(,{t) depends so weakly on t that we can 
neglect the second term in Eq. ()85|) . then in this approximation ( the 
so-called adiabatic approximation) the solution of the evolution equation is 
obviously given by 

a'(t) = e-^^'o^'"(*^'^*a'(to), (86) 

where to is the initial time. 

From Eqs. (|84p and (|86|) for the solution of the evolution equation in the 
flavour representation we have 

a{t) = U'^it) e-'^4 ^'"W'^* f/'"+(to) a{to) . (87) 

Hence, in the adiabatic approximation, the amplitude of the transition Ua — »■ 
Ua' during a time interval (t — to) is given by the expression 

Aiu^ - = E f^r(to) . (88) 

i 

The latter is similar in form to the expression (j39p for the amplitude of 
transitions in vacuum: this is connected with the fact that in the adiabatic ap- 
proximation neutrino remains on the same energy level. The expression (jHHjl 
has a simple meaning: f/^*(to) is the amplitude of the transition from the 
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state of the initial z/q, to the state with energy Eiito) ; the factor e~*''*o ' 
describes the propagation in the state with definite energy; U^-(t) is the am- 
phtude of the transition from the state with energy EY^{t) to the state of the 
final Ua'. The coherent sum over all i must be performed. 

The transition probability must be averaged over the the region where 
neutrinos are produced, over the energy resolution, etc. Oscillatory terms 
in the transition probability usually disappear after the averaging. In this 
case from Eq. (jHHj) we obtain (x ^ t): 



Hence, in the adiabatic approximation the averaged transition probability 
is determined by the elements of the mixing matrix in matter at the initial 
and final points. 

In the case of two neutrino flavours, from Eqs. (f7^ and ()89|) we obtain 
the following expression for the survival probability: 



Let us now consider solar neutrinos in the two-neutrino case. If at some 
point X = xr the MSW [521 ES] resonance condition 



is satisfied, neutrino mixing at this point is maximal {6{xr) = it/4:). 

If the resonance condition is fulfilled, cos26 > (neutrino masses are 
labeled in such a way that Am^ > ). Let us assume that the resonance point 
is in the region between the Sun core, where neutrinos are produced, and the 
surface of the Sun: then at the initial point Xq the electron density will be 
larger than at the resonance point xr and we have A{xo) > Am'^cos26. In 
this case from Eq. (fTHj) it follows that cos 26''" (xq) < and Eq. predicts: 




(89) 




(90) 



Am^ cos 29 = A{xr) 



(91) 



If the condition 



A{xq) » Am^ 



(92) 
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is satisfied, then from Eq. (j76|) we liave 

cos 20™(xo) ~ -1. (93) 

Taking into account tliat 6"^{x) = 6 on the surface of the sun, for the i^e 
survival probabihty we obtain the following relation 

P(z/, ^z/e)~^(l-cos20). (94) 

This relation implies that if the vacuum angle 6 is small, the Uf, survival 
probability is close to zero and practically all z/g's are transfered into other 
neutrino states. 

The MSW resonance condition (jHH) was written in units h = c = 1. It 
can be rewritten as follows 

Am^ cos 2^ ~ 0.7- 10"^^p eV^ (95) 

where p is the density of matter in g ■ cm~^ and E is neutrino energy in 
MeV. In the central region of the Sun p ~ 10^ g ■ cm~^; the energy of solar 
neutrinos is ~ 1 MeV. Thus, the resonance condition is satisfied for the 
solar neutrinos if Arri^ ~ 10~^eV^. 

In the general case of non-adiabatic transitions, for the averaged transi- 
tion probability we have 

P{Ua - Ua') = J2 \U:^kix)\' Pk^ W^ixo)]' (96) 
i,k 

where Pki is the probability of transition from the energy level i to the energy 
level k. 

If we limit ourselves to the two-neutrino case, from the conservation of 
the total probabihty we have 

Pn = I-P21; P22 = I-P12; (97) 

taking into account that P21 = P12, from Eqs. (fTij) . and ^7\\ we obtain 
the following general expression for the z/g survival probability ISHl 

P(z/e Z/e) = ^ + - P12) COS 29'^{x) COS 2^"^(xo) . (98) 
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Different approximate expressions for the transition probability P12 exist 
in tlie literature. In the Landau- Zenner approximation, which is based on 
the assumption that transitions occur mainly in the resonance region, we 
have 
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where 



TT 

e ^ 



Am"^ sin^ 26 







2p cos 26 


1 1^ /^^(^^) 



(99) 



(100) 



For an exponential density F = 1 — tan^ 6; for a linear density F = 1. Let 
us notice that the adiabatic approximation is valid if >> 1. In this case 
P12 ^ 0. 



7 Neutrino oscillation data 

7.1 Evidence in favour of oscillations of atmospheric 
neutrinos 

Atmospheric neutrinos are mainly originated from the decays of charged pions 
and consequent decays of muons: 

n ^ fi + u^] ^ ^ e + + Ue, (101) 

Pions are produced in the processes of interaction of the cosmic rays in the 
atmosphere. In the Super-Kamiokande (S-K) experiment jSHj electron and 
muon neutrinos were detected via the observation of the Cherenkov light in 
the large water- Cherenkov detector (50 kt of II2 O). 

At relatively small energies (< IGeV ) practically all muons decay in 
the atmosphere and from (jlOip it follows that the ratio of the numbers of 
muon and electron events -R^/e must be equal to 2 (if there are no neutrino 
oscillations). At higher energies the expected ratio Rfj,/e is larger than two. 
It can be predicted, however, with an accuracy better than 5 %. 

The ratio -R^/e, measured in the S-K j^H] and S0UDAN2 j^HI atmospheric 
neutrino experiments is significantly smaller than the ratio [Rfi/ e)MC, pre- 
dicted under the assumption of no neutrino oscillations. In the S-K exper- 
iment for the ratio of ratios in the sub-GeV region (-Evis < 1-33 GeV) and 
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Figure 1: Zenith angle distribution of Super-Kamiokande sub-GeV and 
multi-GeV electron and muon events. The black histogram shows the Monte 
Carlo prediction (under the assumption of no neutrino oscillations). The 
gray histogram was obtained for i/^ — > oscillations with the best-fit values 
Am^ = 2.5 X 10-3 eV^ and sin 2d = 1.0. Figure taken from Ref. |i59j . 

multi-GeV region (-Evis > 1.33 GeV) the values 

{R^/e)meas ^ ^ ^ ^ ^ ^ q^^. {Ri,/e)rneas ^ ^ ^ ^ ^ ^ ^^g 

were obtained, respectively. 

The fact that the ratio {R^/e)meas is significantly smaller than the pre- 
dicted ratio {R^/e)MC was known from the results of the previous atmospheric 
neutrino experiments, Kamiokande jHZ| and 1MB |^. For many years this 
"atmospheric neutrino anomaly" was considered as an indication in favour 
of the disappearance of muon neutrinos due to neutrino oscillations. 
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The compelling evidence in favour of neutrino oscillations has been ob- 
tained recently by the S-K collaboration |23, from the observation of the 
large up-down asymmetry of the atmospheric high energy muon events. In 
the S-K experiment [21] the zenith angle dependence of the numbers of elec- 
tron and muon events was measured. In the absence of neutrino oscillations, 
the number of the electron (muon) events in the multi-GeV region must obey 
the following relation: 



where 6z is the zenith angle. 

The results of the S-K experiment on the measurement of the zenith- 
angle distribution of the atmospheric neutrino events is shown in Fig. As 
it is seen from this figure, for the electron events there is a good agreement 
between the data and relation ()102|) . Instead, for the Multi-GeV muon events 
a significant violation of relation (jl(J2|) was observed. The ratio of the number 
of up-going muons (—1 < cos 6z < —0.2) to the number of down-going muons 
(0.2 < cos^^^ < 1) was found to be 



At high energies the direction of leptons practically coincides with the di- 
rection of neutrinos. The up-going muons are produced by neutrinos which 
travel distances from ~ 500 km to ~ 13000 km and the down-going muons 
are produced by neutrinos which travel distances from ~ 20 km to ~ 500 km. 
Hence the observation of the up-down asymmetry clearly demonstrates the 
dependence of the number of muon neutrinos on the distance which they 
travel from the production point in the atmosphere to the detector. 

The S-K data and the data of the other atmospheric neutrino 
experiments, SOUDAN 2 ^ and MACRO [21], are perfectly described if 
we assume that the two-neutrino z/^ — > u^. oscillations take place. From the 
analysis of the S-K data it was found [221 that at 90 % CL the neutrino 
oscillation parameters Am^^^ and sin^ 29atm are in the range 



Aml^ = 2.5 ■ 10^ sin^ 20,tra = 1-0 (xLn = 163.2/170 d.o.f .) (103) 



For the multi-GeV events the effect of the magnetic field of the Earth can be neglected. 



A^/(cos^,) = A^K-cos^,) (/ = e,/i) 



(102) 




0.54 ±0.04 ±0.01. 



1.6 ■ 10^^ < Am^j„ < 3.9 ■ 10^^ eV^; sin^ 2^atm > 0.92; 
the best-fit values of the parameters being equal to 
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7.2 Indications in favour of neutrino oscillations ob- 
tained in the K2K experiment 

Neutrino oscillations in the atmospheric range of Am^ are investigated in 
the first long baseline accelerator experiment K2K [Hni- In this experiment 
neutrinos, originated mainly from the decay of pions, produced at the 12 
GeV KEK accelerator, are recorded by the S-K detector at a distance of 
about 250 km from the accelerator. The average neutrino energy is ~ 1.3 
GeV. 

Two near detectors at a distance of about 300 m from the beam-dump 
target are used in the K2K experiment: a 1 kt water-Cherenkov detector and 
a fine-grained detector. The total number and spectrum of muon neutrinos, 
observed in the S-K detector, are compared to the total number and spectrum 
calculated from the results of the near detectors under the assumption of 
the absence of neutrino oscillations. For the measurement of the energy of 
neutrinos in the S-K detector, quasi elastic one-ring events + n ^ ji" + p 
are selected. 

The first results of the K2K experiment were recently published jHOI- The 
total number of muon events observed in the S-K detector is equal to 56: it 
must be compared with an expected number of events equal to 80.ll5;4. 
The observed number of the one-ring muon events which was used for the 
calculation of neutrino spectrum is equal to 29, while the expected number 
of the one-ring events is equal to 44. 

Thus, in the long baseline accelerator K2K experiment, indications in 
favour of the disappearance of the accelerator were obtained. From the 
maximum likelihood two-neutrino analysis of the data the following best-fit 
values of the oscillation parameters were found 

sin2 2^K2K = l: Am^aK = 2.8 10"^ eV^. (104) 

These values are in agreement with the values of the oscillation parameters 
found from the analysis of the S-K atmospheric neutrino data (see fjl03|) ). 
The first K2K results were obtained with 4.8 10^^ protons on target (POT). 
It is expected that 10^° POT will be utilized in the experiment. 
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7.3 Evidence in favour of the transitions of solar 
into z/^ and Ur 

The energy of the sun is produced in the reactions of the thermonuclear pp 
and CNO cycles, in which protons and electrons are converted into Helium 
and electron neutrinos: 

4p + 2e~ He + 2u^. 

The reactions, which are most important for the solar neutrino experi- 
ments, are listed in Table I. From this Table one can see that the largest part 
of the solar neutrino flux draw up low energy pp neutrinos. According to the 
SSM BPOO [UT], the medium energy mono-energetic ^Be neutrinos make up 
about 10 % of the total flux, while the high energy neutrinos constitute 
only about 10~^ % of the total flux. However, the ^B decay is a very im- 
portant source of the solar neutrinos: in the S-K and SNO jHSl EH HH] 
experiments, due to the high energy thresholds, practically only neutrinos 
from ^B-decay can be detected. The ^B neutrinos also give the dominant 
contribution to the event rate measured in the Homestake experiment |32j . 

Table I 

Reaction Neutrino energy SSM BPOO fiux 

pp de+Ue < 0.42 MeV 5.95 ■ lO^^cm^^s-^ 

e- +^ Be -> z/e ^Li 0.86 MeV 4.77 ■ 10^ cm'^ s"^ 

^B Be* e+ z/e < 15 MeV 5.05 ■ 10^ cm'^ s"^ 



The main sources of the solar neutrinos. The maximum neutrino energies 
and SSM BPOO [HI] fluxes are also given. 

The event rates measured in all solar neutrino experiments are signifi- 
cantly smaller than the ones predicted by the Standard Solar models. 

In the Homestake experiment [321 solar neutrinos are detected through 
the observation of the Pontecorvo-Davis reaction Ue CI e~ +^'' Ar, in 
the GALLEX-GNO jSSj and SAGE |23] experiments through the reaction 
-|-^^ Ga —>■ e~ -f-^^ Ge and in the S-K experiment solar neutrinos are 
detected via the observation of the process z/ + e — > z/ + e. 

^''According to the SSM BPOO the flux of the high energy hep neutrinos, produced in 
the reaction "^He + p He + e+ + 1^0, is about three orders of magnitude smaUer than 
the flux of the neutrinos. 
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For the ratio R of the observed and the predicted by SSM BPOO [HI] rates 
the following values were obtained: 



R = 0.34 ± 0.03 (Homestake) 

R = 0.58 ± 0.05 (GALLEX - GNO) 

i? = 0.60 ±0.05 (SAGE)) 

i? = 0.465 ± 0.018 (S-K). 

If there is neutrino mixing the original solar z/g's, due to neutrino oscilla- 
tions or matter MSW transitions, are transfered into another type of neutri- 
nos, which can not be detected in the radiochemical Homestake, GALLEX- 
GNO and SAGE experiments. In the S-K experiment mainly z/e are detected: 
the sensitivity of the experiment to z/^ and z/^ is about six times smaller than 
the sensitivity to z/g. Thus, neutrino oscillations or MSW transition in matter 
provide a natural explanation for the observed depletion of the fluxes of 
solar Ue- 

Recently a strong model independent evidence in favour of the transition 
of the solar z/g into z/^ and z/^ was obtained in the SNO experiment [HBll37| l3H]. 
The detector in the SNO experiment is a heavy water Cherenkov detector (1 
kton of D2O). Neutrinos from the Sun are detected via the observation of 
the following three reactions: 



1. CC reaction u^ + d^^ +P + P, (105) 

2. NC reaction u + d—^u + n + p, (106) 

3. ES process z/ + e — > z/ + e (107) 

During 306.4 days of running 1967^^^! CC events, 576.5ttli NC events, 
and 263.6!^25 6 events were recorded in the SNO experiment. The kinetic 
energy threshold for the detection of electrons was equal to 5 MeV, the 
NC threshold to 2.2 MeV. Thus, practically only neutrinos from the decay 
'^B — >^ Be + e^ + z/e are detected in the SNO experiment. The important 
point is that the initial spectrum of ^B neutrinos is known [U^ . 
The total CC event rate can be presented in the form 

<^ =< <d > , (108) 
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where < > is the cross section of the CC process (ll(J5|) . averaged over 
the initial spectrum of neutrinos, and ^'^^ is the flux of on the Earth, 
which is given by the relation 

=<P{u,^u,)>cc (109) 

where is the total initial flux of and < P(z/e — > i^e) >cc is the averaged 
z/g survival probability. 

All flavour neutrinos Ue, and Vr are recorded via the detection of the 
NC process ()106|) . Taking into account the z/g — t'^ — Pr universality of the 
neutral current interaction, for the total NC event rate we have 



i?r=<^.T>e^ (110) 

where < a^^ > is the cross section of the NC process ()106p . averaged over 
the initial spectrum of the neutrinos, and is the total flux of all 
flavour neutrinos on the Earth. We have 

E (111) 

l=e,fi,T 

where 

<^ =< P(z/e ^ (112) 

< P(z/e vi) >NC being the averaged probability of the transition — > vi 
{I = e,fi,T). 

All flavour neutrinos are detected also via the observation of the ES pro- 
cess (ll(J7p . However, the cross section of the NC process u^^r + e — > u^^r + e 
is about six times smaller than the cross section of the CC and NC process 

Ue + e ^ Ue + e. 

The total ES event rate can be presented in the form 

=< a.^e > (113) 

Here < a^^e > is the cross section of the process z/gC — > z/gC, averaged over 
initial spectrum of the neutrinos, 

where { ) is the flux of z/^ (u,,. and Ur) and 

0.154. (115) 
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We have 



$f =< P(Z/, ^ Z/,) >E5 (116) 



where < P(z/e — * t^i) >es is the averaged probabihty of the transition ui. 
In the SNO experiment it was found 



('^'f^)sNO = (2.39l|];^^ (stat.) ± 0.12 (syst.)) ■ lO*' cm-^s"^ , (117) 

This value is in a good agreement with the value of the ES flux, obtained in 
the Super-Kamiokande experiment. 

In the S-K solar neutrino experiment [HSI neutrinos are detected via the 
observation of the ES process. During 1496 days of running a large number 
22400 ± 800 of solar neutrino events with recoil total energy threshold of 5 
MeV were recorded. From the data of the S-K experiment the value 

(<l>f^)s-K = (2.35 ± 0.02 (stat.) ± 0.08 (syst.)) ■ 10*^ cm'^g-i . (118) 

was obtained. 

In the S-K experiment the spectrum of the recoil electrons was measured: 
no sizable distortion was observed with respect to the spectrum, expected 
under the assumption of no neutrino oscillations. 

The spectrum of electrons, produced in the CC process p05p . was mea- 
sured in the SNO experiment [38j. Also in this experiment no distortion of 
the electron spectrum was observed. 

Thus, the data of the S-K and SNO experiments are compatible with 
the assumption that in the high-energy region the probability of solar 
neutrinos to survive is practically a constant: 

P(z/e Ue) ^ const. (119) 

The constant z/g survival probability, Eq. (jll9p . implies 

< P{Ue ^e) >CC-< P{j^e ^ J^e) >NC-< P{j^e ^ J^e) >ES ■ 

Taking into account these relations, from (|109p . (jll2p and (|116|) in the high 
energy region we obtain 

. (120) 

In the SNO experiment for the flux of u^. on the Earth it was found 
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Figure 2: Flux of and v^- vs flux of z/g in the '^B energy region deduced 
from the three neutrino reactions in SNO. The diagonal bands show the total 
®B flux as predicted by the BPOO SSM [Hlj (dashed lines) and that measured 
in SNO experiment via the observation of NC (solid band). The interception 
of these bands with the axes represent the ±lcr errors. The bands intersect 
at the fit values for (^^^ and Figure taken from Ref. [HTj . 



($^f )sNO = (1.76+°;°^(stat.)^n9 (syst.)) " lO' cm'^s'^ (121) 

For the flux of all flavour neutrinos (^^^ from the NC measurement it was 
obtained the value 

($DsNO = (5.09l[J:l^(stat.)+S:S (syst.)) ■ 10*^ cm-^s'^ , (122) 

which is about three times larger than the flux of electron neutrinos on the 
Earth. 

Obviously the NC flux is given by 

^ ^NC + ^^c ^123) 

where is the flux of z/^ and is the flux of z/,, and z/^-. 
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By combining now CC and NC fluxes and using the relation (|12(Jj) . we 
can determine the flux <l>^^^. In Refs. jSHl EZl the ES flux (ITTTf) was also 
taken into account (see Fig. For the flux of and on the Earth the 
following value 

(<f..,j5^o = (3.41^0:^^(stat.)^S:S (syst.)) • 10^ cm-^g-i (124) 

was then obtained. Thus, the detection of solar neutrinos via the simultane- 
ous observation of CC, NC and ES processes allowed the SNO collaboration 
to obtain a direct model independent 5.3 a evidence of the presence of 
and Ut- in the flux of the solar neutrinos on the Earth. 

The total flux of the neutrinos, predicted by SSM BPOO jHI], is equal 

to 

(OssMBP = (5.05^^:°}) ■ 10^ cm-^s-i . (125) 

This value is compatible with the value of total flux of all flavour neutrinos 
(|122j) . determined from the data of the SNO experiment. 

The flux of Un and Ur on the Earth can be also obtained from the SNO 
CC data and the S-K ES data. In the first SNO publication |36j the value 

($ .,.Js-K,SNO = (3.69 ± 1.13) • 10^ cm-^s-i , (126) 

was found, which is in a good agreement with the value ()124|) . 

The data of all solar neutrino experiments can be described if one assumes 
that there are transitions of solar into and Ur and that the z/g survival 
probability is given by the two-neutrino expression, which is characterized by 
two oscillation parameters, Am^^^j and tan^ 6'soi. From the global fit of the 
total event rates, measured in all solar neutrino experiments, several allowed 
regions (solutions) in the plane of the oscillation parameters were obtained 
(see, for example, Ref. ^^)'- LMA, LOW, SMA, VO and other regions. 
The situation changed after the recoil electron spectrum was measured in 
the S-K experiment jSni and the SNO results j^ni IHHl were obtained. 
Analyzes of all solar neutrino data showed that the large mixing angle LMA 
MSW region is the most plausible one (see [HI] and references therein). In 
Ref. [23 the following best-fit LMA values of the solar neutrino oscillation 
parameters were found: 

Am^^i = 5 ■ IQ-^eV^; tan^ ^soi = 0.34; xLn = 57/72 d.o.f .. (127) 

Recent data of the KamLAND experiment [311 (see below) allow to ex- 
clude all solutions of the solar neutrino problem, except the LMA one. 
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Figure 3: The ratio of measured to expected numbers of z/e events for different 
reactor neutrino experiments. The shaded region indicates the range of the 
predictions corresponding to the 95% CL. LMA region found from the global 
analysis of the solar neutrino data p^. The dotted curve corresponds to the 
best-fit values Am^^i = 5.5 x 10"^ eV^ and sin^ = 0.83 found in Ref. 
Figure taken from Ref. 

7.4 Reactor experiments CHOOZ and Palo Verde 

In the long baseline reactor experiments CHOOZ jUj and Palo Verde jl2] 
the disappearance of the reactor z/g's in the atmospheric range of Am^ was 
searched for. In spite of the fact that in these experiments no indication in 
favour of neutrino oscillations were found, their results are very important 
for the neutrino mixing. 

In the CHOOZ experiment Ug^s from two reactors at a distance of about 
1 km from the detector were recorded via the observation of the process 

Ue + p ^ + n. 

The value of the ratio R of the total number of the detected events to 
the expected number was found to be 

R = 1.01 ± 2.8% (stat) ± ±2.7% (syst) (CHOOZ) 
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Figure 4: Upper panel: Expected reactor z/g energy spectrum with contribu- 
tions of z/geo (antineutrinos produced in ^ss-q ^^^^ "^^"^Th decays in the Earth) 
and accidental background. Lower panel: Energy spectrum of the observed 
prompt events (solid circles with error bars), along with the expected no os- 
cillation spectrum (upper histogram, with z/geo and accidentals shown) and 
best fit (lower histogram) in the case of neutrino oscillations. The shaded 
band indicates the systematic error in the best-fit spectrum. The vertical 
dashed line corresponds to the threshold of 2.6 MeV. Figure is taken from 
Ref. [SH]. 

In the similar Palo Verde experiment it was obtained: 

R = 1.01 ± 2.4% (stat) ± 5.3% (syst) (Palo Verde) 

The data of the experiments were analyzed in ^T] 02] in the framework 
of two-neutrino oscillations and exclusion plots in the plane of the oscillation 
parameters Am^ and sin^ 2 9 were obtained. From the CHOOZ exclusion 
plot at Am^ = 2.5 ■ 10"^ eV^ (the S-K best-fit value) one gets 

(sin2 2^)cHooz <1.5-10-^ 
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7.5 The KamLAND evidence in favour of the disap- 
pearance of the reactor z/g 

The first results of the KamLAND experiment [221 were published recently. 
In this experiment z/g's from many reactors in Japan and Korea are detected 
via the observation of the classical process 

Ue + p ^ e'^ + n 

The threshold of this process is ~ 1.8 MeV. About 80 % of the total number 
of the events is due to Ue from 26 reactors within the distances of 138-214 
km. 




Figure 5: KamLAND excluded regions of neutrino oscillation parameters 
Am^ and sin^ 29 for the rate analysis and allowed regions for the combined 
rate and energy spectrum analysis at 95% C.L. At the top are the 95% 
CL. excluded region from CHOOZ jJT] and Palo Verde [32] experiments, 
respectively. The dark area is the 95% C.L. LMA allowed region obtained 
in Ref. |65^. The thick dot indicates the best-fit values of the oscillation 
parameters. Figure taken from Ref. |39j . 

The 1 kt liquid scintillator detector of the KamLAND experiment is lo- 
cated in the Kamioka mine at a depth of about 1 km. Both prompt photons 
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from the annihilation of e"*" in the scintillator and 2.2 MeV delayed photons 
from the neutron capture n + p ^ c? + 7 are detected. The mean neutron 
capture time is 188 ± 23 /isec. In order to avoid background, mainly from 
the the decay of ^^^U and ^^^Th in the Earth, the cut Eprompt > 2.6 MeV was 
applied. 

During 145.1 days of running there were observed 54 z/g events. The 
number of the events expected in the absence of neutrino oscillations is 
equal to 86.8 ± 5.6. For the ratio of observed and expected Ue events the 
following value 

^^^-^^^ = 0.611 ± 0.085 ± 0.041. (128) 

was obtained. 

In Fig. ini for all reactor neutrino experiments, the dependence of the 
ratio of the observed and expected events on the average distance between 
reactors and detectors is plotted. The dotted curve was calculated with 
the best-fit solar neutrino LMA values of the oscillation parameters Am^ = 
5.5 ■ lO-^eV^ and sin^ 26 = 0.83, obtained in Ref. [03 ■ 

In the KamLAND experiment the prompt energy spectrum was also mea- 
sured (see Fig. E}. The prompt energy is connected with the energy of Ue 
by the relation Eprompt = Ep^ — 0.8 MeV — En {E^ being the average energy 
of the neutron). From the two-neutrino analysis of the KamLAND data the 
following best-fit values of the oscillation parameters were obtained 

(Am^)KamLAND = 6.9 ■ 10~^eV^; (sin^ 2 6')KamLAND = 1 • 

These values are compatible with the values of the oscillation parameters 
in the solar neutrino LMA MSW region. The allowed region in the plane of 
the oscillation parameters, obtained from the analysis of the measured rate 
and measured spectrum is shown in Fig. El (95 % CL). The region outside 
the solid line is excluded from the rate analysis. The dark region is the solar 
neutrino LMA allowed region, obtained in Ref. 

The KamLAND results provide a strong evidence for neutrino masses 
and oscillations, obtained for the first time in an experiment with terrestrial 
antineutrinos with the expected flux well under control. It allows us to 
exclude the SMA, LOW and VAC regions of neutrino oscillation parameters. 
The only viable solution of the solar neutrino problem appears to be the 
LMA MSW solution. 
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8 Neutrino oscillations in the framework of 
three-neutrino mixing 

8.1 Neutrino oscillations in the atmospheric range of 
the neutrino mass-squared difference 

We have discussed evidences in favour of neutrino oscillations, which were 
obtained in the solar, atmospheric and reactor KamLAND neutrino experi- 
ments. At present there exists also an indication in favour of the transition 
i^/i ^e, that was obtained in the single accelerator LSND experiment 
The LSND data can be explained by neutrino oscillations. From the analysis 
of the data, the following ranges for the values of the oscillation parameters 



were obtained. 

In order to describe the data of the solar, atmospheric, KamLAND and 
LSND experiments, which require three different values of neutrino mass- 
squared differences, it is necessary to assume that there are (at least) four 
massive and mixed neutrinos. This means that in addition to the three 
flavour neutrinos (at least) one sterile neutrino must exist (see, for example, 
Ref. i35|)- The result of the LSND experiment requires, however, confirma- 
tion. The MiniBooNE experiment at Fermilab [HZI, which started in 2002, 
is aimed at checking the LSND result. 

We will consider here the minimal scheme of three-neutrino mixing 



where U is the unitary 3x3 PMNS mixing matrix. This scheme provides two 
independent Am^'s and allows us to describe solar, atmospheric, KamLAND 
and other neutrino oscillation data. 

Let us start with the consideration of neutrino oscillations in the atmo- 
spheric range of Am^, which can be explored in the atmospheric and long 
baseline accelerator and reactor neutrino experiments. In the framework of 
the three-neutrino mixing ()129|) . with mi < m2 < rris, there are two possi- 
bilities: 



2- 10-i<(Am2)LSND<leV2; 

3- 10-3<(sin2 2^)LSND<4-10-2 



3 




(129) 



1=1 
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I. Hierarchy of neutrino mass- squared differences 

Am^i ~ Ami,- Am^^ ^ Aml^; Am^, < Aml^. (130) 

II. Inverted hierarchy of neutrino mass-squared differences 

Am^^ ~ Am^^i; Am2, ^ Am^,^; Am^^ < Am^,. (131) 

We will first assume that the neutrino mass spectrum is of the type I. 
The values of L/E, relevant for neutrino oscillations in the atmospheric range 
of neutrino mass-squared difference, satisfy the inequality: 



Am2i ^ < 1. 



L 
E 

Thus, we can neglect the contribution of Am2i to the transition probability, 
Eq. ()40|) . In this case, for the probability of the transition i^a' we 

obtain the following expression 

2 



Sa'a + U^'sK, (e-^^'^i^^ - 1) . (132) 



Hence, in the leading approximation, the transition probabilities in the 
atmospheric range of Am^ are determined by the largest neutrino mass- 
squared difference Am'^2 ^"^^ by the elements of the third column of the 
neutrino mixing matrix, which connect the flavour neutrino fields with 
the field of the heaviest neutrino u^l. 

For the appearance probability we obtain from Eq. (jl32j) the expression 



P{Ua l^a') = ^^a';a (^^ - COS Ar/iga (a 7^ a'), (133) 

where the oscillation amplitude is given by 

A„';„ = 4 |t/,,3|' |f/a3|'. (134) 

The survival probability can be obtained from the condition of conserva- 
tion of probability and from Eqs. p33j) and p 331 we have 



p(z/« i^«) = 1 - 5Z ~^ ^""'^ y- ~ cos Am32 ^ j 



(135) 
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Taking into account the unitarity of the mixing matrix, the amphtude Bq-^ 
can be written as: 

Ba;« = = 4 |f/„3|' (1 - \Ua3\^) . (136) 

Obviously, in the case of inverted hierarchy of the neutrino mass-squared 
differences (case II above) the transition probabilities can be obtained from 
Eqs. fll33|) - p36|) with the replacements Am^g ^^21 If/asp — » |f/aip. 

We notice that the transition probability ()133p depends only on It/asp and 
Am32. The CP phase does not enter into this equation and in the leading 
approximation the relation 

P(Z/^ Ua') = P{Ua i^a') (137) 

is automatically satisfied. 

Thus, the investigation of the effects of CP violation in the lepton 
sector in the future long baseline neutrino oscillation experiments will be a 
difficult problem: possible effects are suppressed due to the smallness of the 
parameter Am^2/^"^32- High precision experiments on the search for effects 
of CP- violation in the lepton sector are planned for the future neutrino facility 
JHF [68J and Neutrino Factories 

The transition probabilities ()133p and p35|) have a two-neutrino form 
in each channel. This is the obvious consequence of the fact that only the 
largest mass-squared difference Am^g contributes to the transition proba- 
bilities. The elements |t/a3p, which determine the oscillation amplitudes, 
satisfy the unitarity condition |^a3p = 1- Hence, in leading approxima- 
tion, transition probabilities are characterized by three parameters. We can 
choose the latter to be: 

Am32, sin^6'23, |f/e3p 

Then, from Eqs. and ()134|) . for the amplitudes of z/^ — >■ and z/^ — >■ i^e 
transitions we obtain the expressions 

^.;^ = (l-|t/e3r)'sin2 2^23; 

^;/. = 4|[/e3r(l-|f/e3|'sin2 ^23- 
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8.2 Oscillations in the solar range of neutrino mass- 
squared difference 

Let us consider now, in the framework of the tree-neutrino mixing, neutrino 
oscillations in the solar range of Am^. The z/g survival probability in vacuum 
can be written in the form 



-i Amf 



^ + \Ue3\^e 



-i Am: 



31 2E 



138) 



«=1,2 



We are interested in the z/g survival probability averaged over the region where 
neutrinos are produced, over energy resolution etc. Because of the hierarchy 
Am32 3> Am2i, in the expression for the averaged survival probability the 
interference between the first and the second terms in Eq. (jl38|l disappears. 
The averaged survival probability can then be presented in the form 



Here P^^'^^u^ 
where 



|t/e3r+(l-|f/e3p)' P^^'^^ (^e - ^e) 



z/g) is given by the expression 
1 



1 



^(1,2) 



^ 2 L\ 
1 — COS Am 



''21 



2E J 



\Uel\ 


to 


Ue2 


2 


(1- 


\Ue3 


2 


to 1 



With the help of Eq. (gSl), Eq. (|TTr]l simplifies to 

A(i'2)=sin2 2^i2. 



(139) 

(140) 
(141) 

(142) 



Thus, the probability P*^^'^)(z/e — > z/g) is again characterized by two parame- 
ters only and has the standard two-neutrino form. 

The expression fll39|) is also valid for oscillations in matter (see j7T| I45j). 
In this case P'^^'^^(z^g — > z^g) is the two-neutrino z/g survival probability in 
matter, calculated under the condition that the density of electrons pg(x) 
in the effective Hamiltonian of the interaction of neutrino with matter is 
changed by (1 - |t/e3p) pg(x). 
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Hence, the Ve survival probability is characterized, in the solar range of 
Am^, by three parameters 



Am2i, tan^6'i2, |f/, 



|2 



e3| 

The only common parameter for the atmospheric and solar ranges of Am^ 
is |f/e3p- As we will see in the next subsection, from the data of the reactor 
CHOOZ and Palo Verde experiments, this parameter turns out to be small. 

8.3 The upper bound of \Uei\'^ from the CHOOZ data. 
DecoupUng of oscillations in the solar and atmo- 
spheric ranges of Am^ 

The long baseline reactor experiments CHOOZ |1T] and Palo Verde are 
sensitive to the atmospheric range of Am^. No indications in favour of the 
disappearance of reactor z/g was obtained in these experiments. From the 
analysis of their data a stringent bound on the parameter |?7e3p was ob- 
tained, as it is explained below. 

In the framework of the three-neutrino mixing the probability of z/g to 
survive in the atmospheric range of Am? is given by the expression 

P(z7e ^ Z7,) = 1 - iPe^e {l " COS Am^^ , (I43) 

where 

Be;e = 4|?7,3p (l-|f/e3|') • (144) 

In Refs. j4H W]\ exclusion plots in the plane of the parameters Am^ = Am\2 
and sin^ 26 = Bg-e were obtained, from which one has 

Be;e < B°„ (145) 

where the upper bound B^.g depends on the S-K allowed values 

of the parameter Amgg, from the CHOOZ exclusion plot we have 

1 ■ 10^^ < B°.g < 2.4 ■ 10"^ (146) 



From Eqs. ()144|) and ()146p the following bounds on the parameter \Ues\ 
can be easily obtained: 

|[/g3r < - a/1 - Ke) ^ 7 Ke (147) 



2 \ V e;e y ~ 4 
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or 

Thus, the parameter |t/e3p can be either small or large (close to one). 
This last possibility is excluded, however, by the solar neutrino data. In 
fact, if |t/e3p is large, from Eq. it follows that in the whole range of 

the solar neutrino energies the probability of z/g to survive is close to one in 
obvious contradiction with the solar neutrino data. Hence, the upper bound 
of the parameter |?7e3p is given by Eq. ()147|) . At the S-K best-fit value 
Am|2 = 2.5 ■ lO-^eV^ we get: 

If/esl^ < 4 ■ 10-2 (95% CL). (149) 

Taking into account the accuracies of the present-day experiments, one 
can neglect |?7e3p in the expressions for the transition probabilities. In this 
approximation, neutrino oscillations in the atmospheric range of Am^ are 
oscillations, with Am^j — ^^Itm Ar-fj, — sin^ 2 6*23 — sin^ 2 6'atm- 
In the solar range of Am^ we have 

P(z/e^Z/e)~P(l'2)(z/, -.z/e), (150) 

where the two-neutrino survival probability P(z/e — * z/g) depends on the 
parameters Amgi — Am^^j and tan^ 2 612 ~ tan^ 2 ^^soi- 

Thus, due to the smallness of the parameter |f/e3p and the hierarchy 
of neutrino mass squared differences Am^2 <S Arrigg neutrino oscillations in 
the atmospheric and solar ranges of Am^, in the leading approximation, are 
decoupled [721 ES] and are described by the two-neutrino formulas, which 
are characterized, respectively, by the oscillation parameters 

Am32, sin2 2^23 and Am2i, tan2^i2 

. From the CHOOZ data only the upper bound of the parameter |?7e3p 
can be obtained. The possibilities to investigate effects of the three-neutrino 
mixing and in particular important effects of CP violation in the lepton sector 
depend on the value of this parameter. The value of the parameter |t/e3p 
will be probed in MINOS [73] and ICARUS [73] experiments and in neutrino 
experiments at JHF [73] and Neutrino Factories (see [U^ I70j ). 

From the data of the experiments on the investigation of neutrino os- 
cillations only neutrino mass-squared differences Amgi and Amgg can be 
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determined. Neutrino masses m2 and rris are given by the relations 



777,2 = Y + Am^i ; ms = y + Amgg + Am2i . (151) 

In the next section we will discuss the resuhs of experiments which allow 
us to obtain information about the absolute values of neutrino masses. 



9 p-decay experiments on the measurement 
of neutrino mass 

The standard method for the measurement of the absolute value of the 
neutrino mass is based on the detailed investigation of the high-energy part 
of the /3-spectrum of the decay of tritium 

^3 He + e- + Ue. (152) 

This decay is the super-allowed one. Thus, the nuclear matrix element is a 
constant and the electron spectrum is determined by the phase space. The 
decay (jl52p has a small energy release {Eq ^ 18.6 keV) and a convenient 
time of life (T1/2 ~ 12.3 years). 

The standard effective Hamiltonian of the /3-decay is given by 

nf'' = ^ 2ei7„z/,i r + h.c. . (153) 

Here is the hadronic charged current, Gp is the Fermi constant and 

^eL = ^ Ueil^iL , (154) 
i 

where i/j is the field of neutrino with mass rrii and U is the unitary mixing 
matrix. 

Neglecting the recoil of the final nucleus, from Eqs. p53|l and ()154|1 for 
the spectrum of electrons we obtain the following expression 

^ = ym'^, (155) 
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where 



dTi 



(156) 



dE 



Cp{E + me) {Eo - E) ^{Eo-EY-mfFiE) d{E^ - E - rm) , 



E being the kinetic energy of the electron, Eq the energy released in the 
decay and rrie is the mass of the electron. The Fermi function F{E) takes 
into account the Coulomb interaction of the final particles and the constant 
C is given by the expression 



where ^ c is the Cabibbo angle and M the nuclear matrix element ( which 
is a constant). 

The sensitivity to the neutrino mass of the present-day tritium experi- 
ments Troitsk [ZE] and Mainz j2Z| is 2-3 eV. In the future experiment KA- 
TRIN [ZH! a sensitivity of 0.25 eV is expected. Taking into account that 
-^/Amf^ ~ 5 lO^^eV is much smaller than the sensitivities of the present 
and future tritium experiments, from the relations we can conclude 

that the neutrino mass can be measured in these experiments only if the 
neutrino mass spectrum is practically degenerate: mi ~ m2 — m^. In this 
case, from the unitarity of the mixing matrix we have: 



Let us discuss now the results of the tritium experiments. In the Mainz 
experiment [ZZj the target is molecular tritium condensed on a graphite sub- 
strate. The spectrum of the electron is measured by the integral electrostatic 
spectrometer, which combines high luminosity with high resolution. The res- 
olution of the Mainz spectrometer is equal to 4.8 eV. In the analysis of the 
experimental data four free variable parameters are used: the normalization 
C, the background B, the released energy Eq and the neutrino mass-squared 
ml. The analysis of 1998, 1999 and 2001 data gave the following result: 



C = Gj, 




cos^ ec\M\'^ 




(157) 



mj = (-1.2 ±2.2 ±2.1) eV^ 



(158) 



This value corresponds to the upper bound 



mi < 2.2 eV (95% CL) . 



(159) 
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The integral electrostatic spectrometer is also used in the Troitsk neutrino 
experiment JTB]: its resolution is 3.5-4 eV. In the Troitsk experiment the 
source is a gaseous molecular tritium source. From the four-parameter fit 
of the Troitsk data, for the parameter mf large negative values, in the range 
(—10 -. — 20) eV^, were obtained. The investigation of the character of the 
measured spectrum suggests that the effect of obtaining a negative ml is 
due to a step function superimposed on the integral continuous spectrum: 
the step function in the integral spectrum corresponds to a narrow peak in 
the differential spectrum. 

In the analysis of the data, the authors of the Troitsk experiment added to 
the theoretical integral spectrum a step function with two additional variable 
parameters ( the position of the step £^step and its height). From the six- 
parameter fit of the data, the following value for the parameter ml 

m\ = (-2.3 ± 2.5 ± 2.0) eV^ . (160) 

was found. From ()160|) the upper bound on the neutrino mass mi 

mi < 2.2 eV (95% CL) . (161) 

can be deduced. 

The analysis of the Troitsk data shows that the position of the step Eq — 
Estep is periodically changed in the interval 5-15 eV and the average value of 
the height of the step is about 6 • 10^^^. The existence of this anomaly was 
not confirmed by the Mainz experiment [77]. 

A new tritium experiment, KATRIN [7H], is now under preparation. 
In this experiment gaseous molecular source and frozen tritium source are 
planned to be used. The integral electrostatic spectrometer will have two 
parts: the pre-spectrometer, which will select electrons in the last ~ 100 eV 
of the spectrum, and the main spectrometer. The latter will have a res- 
olution of ~ 1 eV. It is expected that the KATRIN experiment will start 
to collect data in 2007. After three years of running the accuracy of the 
measurement of the neutrino mass will reach ~ 0.25 eV. 

In order to understand the origin of the small neutrino masses we need 
to know the nature of the massive neutrinos: are they Majorana or Dirac 
particles? As we have seen, neutrino oscillation experiments can not answer 
this fundamental question. The nature of the massive neutrinos can be re- 
vealed in experiments on the search for neutrinoless double /3-decay of some 
even-even nuclei. The next section will be devoted to the discussion of such 
experiments. 
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10 Neutrinoless double /3-decay 



The search for neutrinoless double /?-decay 

{A,Z) ^ {A,Z + 2) + e- + e- (162) 

of some even-even nuclei is the most sensitive and direct way of investigating 
the nature of neutrinos with definite masses. The total lepton number in the 
process (I162j) is violated and this is allowed only if the massive neutrinos 
Vi are Majorana particles. 

We will assume that the Hamiltonian of the process has the standard 
form, Eq. ()153j) . and that the flavour field PeL is given by 

T^eL = Uei^iL , (163) 

i 

where Ui are Majorana fields. 

The neutrinoless double /3-decay {{P P)ou -decay) is a process of second 
order in the Fermi constant Gp, with virtual neutrinos. For small neutrino 
masses the neutrino propagator is given by the expression 

<0\T{u^Uxi)^J,{x,))\0>^m,, ^ J ^e-^''(---)i^C, (164) 

where 

^PP = '^U^i^i- (165) 

i 

The total matrix element of the {PP)ou -decay is a product of and 
the nuclear matrix element, which does not depend on neutrino masses and 
mixing. 

Results of many experiments on the search for (/5/3)o,y-decay are available 
at present (see Refs. [TnHHOl)- No indication in favour of (/?/3)oy-decay was 
obtained so far.^^ The most stringent lower bounds for the time of life of 
(/3/3)oi/-decay were obtained in the Heidelberg-Moscow [84j and IGEX jH3] 
''^Ge experiments: 

T°;2 > 1.9 ■ 10^^ y (90% CL) Heidelberg - Moscow 

> 1-57 ■ 10^^ y (90% CL) IGEX 

^'^Thc recent claim |H3 of evidence of a {(3/3)o ,y-decay, obtained from the reanalysis of the 
data of the Heidelberg-Moscow experiment, has been strongly criticized in Refs. |82l ISHj . 



51 



Taking into account different calculations of the nuclear matrix element, from 
these results the following upper bounds were obtained for the effective Ma- 
jorana mass: 

\mpp\ < (0.35- 1.24) eV. (166) 

Many new experiments on the search for the neutrinoless double /3-decay are 
in preparation at present (see Ref. [70]). In these experiments the sensitivi- 
ties 

\mpp\ ^ (1 ■ 10"^ - 1.5 ■ 10-2) eV 

are expected to be achieved. 

The evidence for neutrinoless double (3- decay would be a proof that neu- 
trinos with definite masses are Majorana particles. The value of the effective 
Majorana mass \'mjip\ combined with the values of the neutrino oscillation 
parameters, obtained from the results of neutrino oscillation experiments, 
would enable us to obtain important information about the character of the 
neutrino mass spectrum, the minimal neutrino mass mi and the Majorana 
CP phase (see Ref. jHE] and references therein). 

We will consider here three typical neutrino mass spectra. 

1. The hierarchy of neutrino masses mi <^ m2 <^ m^. 

In this case for the effective Majorana mass we have the following upper 
bound : 

\m^p\ < sin^ ^,,1 ^Am^^i + \Uez? ^A^^atm- (167) 

Using the best-fit values of the oscillation parameters and the CHOOZ 
limit on |f/e3p (see Eqs. (jS2I), dZZD and (fT10|l ). we obtain from (jTTTTj) 
the upper bound 

I < m > I < 3.8 ■ 10"^, (168) 

which is significantly smaller than the expected sensitivities of the fu- 
ture {PP)ou- experiments. 

2. Inverted hierarchy of neutrino masses: mi <^ m2 < m^. 

The effective Majorana mass is given, in this case, by the expression 

|"i/3/3| ^ (l - sin2 2 6'soi sin^a)^ ^y Aml^^, (169) 

where a = — a2 is the the difference of the Majorana CP phases 
{Uei = \Uei\ e*°'). From this expression it follows that 

V^Amf^ I cos 2e,,i\<\mpp\<^Aml^, (170) 
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where the upper and lower bounds correspond to the case of CP conser- 
vation with equal and opposite CP parities of 1/3 and z/2, respectively. 

Using the best-fit value of the parameter taiP 9soi (see Eq. (|127p ). we 
have 

1 



2 , Aml^ < \mp^\ < V^rnl,^ , (171) 

Thus, in the case of the inverted mass hierarchy, the scale of \mpi3\ is 
determined by A/Amf^. Should the value of be in the range 

(jlTlj) . (which can be reached in the future experiments on the search 
for (/5/3)o!/-decay), then we would have an argument in favour of the 
inverted neutrino mass hierarchy. 

3. Practically degenerate neutrino mass spectrum: m2 — ~ mi ^ 

The effective Majorana mass in this case is given by the expression 

3 

K^l^mi \J2U!.\- (172) 

i=l 

Neglecting the small contribution of |?7e3p ( l^eiP in the case of the 
inverted hierarchy), for one obtains relations Eqs. ()169p - ()171|) in 

which a/AtttI^ is replaced by mi. Thus, a signature of the degenerate 
neutrino mass spectrum would be the occurrence of a value |m^^| ^ 

It is obvious that for neutrino mass mi we have the following bound 

\mpp\ < mi < '""fj ^ 2 Im^^l (173) 

I cos / Usol I 

and the parameter sin^ a, which characterizes the violation of the CP 
invariance in the lepton sector (HH], is given by the relation 

Thus, if the mass mi is measured in the future experiments and the 
value of the parameter sin^ 26'soi is determined with high accuracy in 
the KamLAND BOREXINO [HHl and other neutrino experiments, 
from the results of the future (/9/9)oi/ experiments information on the 
Majorana CP phase can be inferred. 
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All previous conclusions are based on the assumption that the value of 
the effective Majorana mass can be obtained from the measurement of 

the life-time of the (/3/3)oi/-decay. However, the determination of the param- 
eter \Tnpp\ from the experimental data requires the knowledge of the nuclear 
matrix elements. At present there are large uncertainties in the calculation of 
these quantities (see, for example, Refs. [HSl EOl EI] ) ■ Different calculations 
of the lifetime of the (/?/3)oj.-decay differ by about one order of magnitude. 

In Ref. a method was proposed, which allows one to check the results 
of the calculations of the (/3/3)oj/-decay nuclear matrix elements in a model 
independent way. This method can be applied if (/3/3)o;y-decay of different 
nuclei is observed. 

11 Conclusion 

About forty years after the original idea of B. Pontecorvo, compelling evi- 
dence in favour of neutrino oscillations were obtained in the S-K [29j, SNO fHHl 
EHl EZl; KamLAND [HH] and other neutrino experiments. These findings 
opened a new field of research in particle physics and astrophysics: the physics 
of massive and mixed neutrinos. 

From the results of the experiments it follows that neutrino masses are 
many orders of magnitude smaller than the masses of leptons and quarks. 
Tiny neutrino masses are a first signature of a new physics, beyond the 
Standard Model. 

There remain many unsolved problems in the physics of massive and 
mixed neutrinos. The problem concerning |?7e3p is an urgent one: on 
the value of |?7e3p depends the possibility to study the effects of the three- 
neutrino mixing and in particular the effects of CP violation in the lepton 
sector. 

Another problem is the one connected with the LSND experiment |66j . 
If the LSND result will be confirmed by the MiniBOONE experiment |U7] . 
this will imply that the number of light neutrinos is more than three and in 
addition to the three flavour neutrinos sterile neutrino (s) must exist. If, on 
the contrary, LSND result is refuted, the minimal scheme with three massive 
and mixed neutrinos will be a plausible possibility. 

The problem of the nature of massive neutrinos (Dirac or Majorana?) is 
the most fundamental one. This problem can be solved by the experiments on 
the search for neutrinoless double /3-decay. From the existing data the foUow- 
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ing bound on the effective Majorana mass was found: Impjjl < (0.3 — 1.3) eV. 
In future experiments, now in preparation, a significant improvement of tlie 
sensitivity is expected. 

Due to tlie interference nature of tfie neutrino oscillation phenomenon 
and to the possibilities of exploring large values of L/E, neutrino oscillation 
experiments are sensitive to very small values of the neutrino mass-squared 
differences. The determination of the absolute value of the neutrino masses 
requires the so-called direct measurements and it is a challenging problem. 
From the data of the tritium experiments the upper bound mi < 2.2 MeV 
was obtained. The future experiment KATRIN [7S] is expected to be sensitive 
to mi ~ 0.25 MeV. 

The progress in cosmological measurements, achieved in the last years, 
allows one to reach ~ 1 eV sensitivity for the sum of the neutrino masses 
The most stringent bound. 



was recently obtained in Ref. jHl]. For a degenerate neutrino mass spectrum 
this bound implies 



In Ref. |91| the 2dF Galaxy Redshift Survey ^ data, recent WMAP high 
precision data and other cosmological data were used. A significant improve- 
ment of this limit is expected with the future Sloan Digital Sky Survey jHS] , 
future WMAP and other data. 

We have discussed here mainly the phenomenology of neutrino mixing 
and oscillations, and the most recent experimental data. For models of neu- 
trino masses and mixing we refer the reader to the recent reviews |46| ITf| HH] . 
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